Citrus-greening (huanglongbing)-induced small rnas are potential early diagnosis markers

ABSTRACT

The present invention provides compositions and methods for detecting  Candidatus Liberibacter  infection and Huanglongbing disease in a citrus plant by detecting the expression of small RNAs such as miRNA and siRNA. The invention also provides methods for treating Huanglongbing disease in a citrus plant by contacting the plant with a phosphorus containing solution.

CROSS-REFERENCES TO RELATED APPLICATIONS

The present application is a Continuation of U.S. patent applicationSer. No. 14/003,760, filed Dec. 9, 2013, which claims benefit ofpriority to International Application No. PCT/US2012/028104, filed Mar.7, 2012, and U.S. Provisional Patent Application No. 61/450,050, filedMar. 7, 2011, the contents of each of which are incorporated herein byreference in their entirety.

REFERENCE TO A “SEQUENCE LISTING,” A TABLE, OR A COMPUTER PROGRAMLISTING APPENDIX SUBMITTED AS AN ASCII TEXT FILE

The Sequence Listing written in file 081906-1057903-211120US.txt,created on Aug. 4, 2017, 22,288 bytes, machine format IBM-PC, MS-Windowsoperating system, is hereby incorporated by reference in its entiretyfor all purposes.

BACKGROUND OF THE INVENTION

Citrus greening, also called Huanglongbing or yellow dragon disease, isa disease of citrus. This bacterial disease is thought to haveoriginated in China in the early 1900's. The disease is primarily spreadby two species of psyllid insects. One species, the Asian citruspsyllid, Diaphorina citri, has been present in Florida since 1998. Thebacteria itself is not harmful to humans but the disease is damaging tothe citrus crops. There are three strains of the bacteria: an Asianversion, an African version, and a recently described American straindiscovered in Brazil.

The Asian strain, Candidatus Liberibacter asiaticus (Ca. L., Las), wasfound in Florida in early September, 2005. As a result, citrus greeningdisease is becoming a major threat to the U.S. citrus industry. Otherthan tree removal, there has been no known effective control once a treeis infected and there has been no known cure for the disease. Infectedtrees may produce misshapen, unmarketable, bitter fruit, or no fruit.Citrus greening reduces the quantity and quality of citrus fruits,eventually rendering infected trees useless. In areas of the worldaffected by citrus greening the average productive lifespan of citrustrees has dropped from 50 or more years to 15 or less. The trees in theorchards usually die 3-5 years after becoming infected and requireremoval and replanting. An infected tree produces fruit that isunsuitable for sale as fresh fruit or for juice.

Citrus plants infected by the citrus greening bacteria may not showsymptoms for years following infection. Initial symptoms frequentlyinclude the appearance of yellow shoots on a tree. As the bacteria growand move within the tree, the entire canopy progressively develops ayellow color.

The most characteristic symptoms of citrus greening are a blotchy leafmottle and vein yellowing that develop on leaves attached to shoots,providing the overall yellow appearance. These foliar symptoms maysuperficially resemble a zinc deficiency although the green and yellowcontrast is not as vivid with greening as it is with zinc deficiency oranother disease, citrus variegated chlorosis. Leaves with citrusgreening have a mottled appearance that differs from nutrition-relatedmottling in that greening-induced mottling usually crosses leaf veins.Nutrition related mottles usually are found between or along leaf veinsand leaves may be small and upright.

Fruit from diseased trees are small, often misshapen, and typically somegreen color remains on ripened fruit. On Mandarin orange, fruit maydevelop an uneven ripening such that they appear half orange and halfyellow. This symptom is the origin of the common name “greening.” Yieldsare almost minimal, and any developed fruit is rendered worthless due tosmall size, poor color, and bad taste.

BRIEF SUMMARY OF THE INVENTION

The present disclosure provides compositions and methods for determiningif a citrus plant has Huanglongbing (HLB) disease or is infected withCandidatus Liberibacter (Ca. L.). In one aspect, HLB disease orinfection with Ca. L. is detected by detecting the expression level of asmall RNA, such as a micro RNA or small interfering RNA, in the plant,and comparing the level of the small RNA to the level detected in acontrol, non-infected, plant. Thus, in some embodiments, a method fordetection of HLB disease or Ca. L. asiaticus-infection in a citrus plantis provided, the method comprising detecting in a sample from the citrusplant the level of expression of one or more RNAs selected from thegroup consisting of miRNA399, siRNA1005, siRNA1008 or siRNA1009, whereinincreased expression of the one or more RNAs compared to expression ofthe RNA in a non-infected healthy plant indicates the presence of Ca. L.asiaticus-infection.

In some embodiments, miRNA399 comprises the consensus sequence TG X₁ X₂AAAGGAG X₃ X₄ TTGCC X₅ X₆ X₇, where X₁ is C or T, X₂ is C or T, X₃ is Aor C, X₄ is G, T, or A, X₅ is C or A, X₆ is T or G, and X₇ is A or G(SEQ ID NO:99).

In some embodiments, miRNA399 comprises a sequence that is substantiallyidentical to (e.g., at least 80%, 85%, 90%, or 95% identical to) asequence selected from TGCCAAAGGAGAGTTGCCCTG (SEQ ID NO:1),TGCCAAAGGAGAGTTGCCCTA (SEQ ID NO:2), TGCCAAAGGAGATTTGCCCGG (SEQ IDNO:8), and TGCCAAAGGAGAATTGCCCTG (SEQ ID NO:6).

In some embodiments, miRNA399 comprises a sequence that is substantiallyidentical to (e.g., at least 80%, 85%, 90%, or 95% identical to) asequence selected from UGCCAAAGGAGAUUUGCCCGG (SEQ ID NO:9),UGCCAAAGGAGAGUUGCCCUA (SEQ ID NO:10), UGCCAAAGGAGAAUUGCCCUG (SEQ IDNO:11), and UGCCAAAGGAGAGUUGCCCUG (SEQ ID NO:12).

In some embodiments, the siRNA1005 comprises a sequence that issubstantially identical to (e.g., at least 80%, 85%, 90%, or 95%identical to) the sequence ATAGATAATGGATCAACGGTTATA (SEQ ID NO:13). Insome embodiments, the siRNA1008 comprises a sequence that issubstantially identical to (e.g., at least 80%, 85%, 90%, or 95%identical to) the sequence TCGAACAAGGTAAGGATGTCA (SEQ ID NO:14). In someembodiments, the siRNA1009 comprises a sequence that is substantiallyidentical to (e.g., at least 80%, 85%, 90%, or 95% identical to) thesequence CTTCTAATAAACATGCATGAA (SEQ ID NO:15).

In some embodiments, the method further comprises detecting the mRNA ofa ubiquitin-conjugating enzyme gene. In one embodiment, theubiquitin-conjugating enzyme is ubiquitin-conjugating enzyme E2 (UBC).In one embodiment, the UBC gene or mRNA comprises a sequence that issubstantially identical to (e.g., at least 80%, 85%, 90%, or 95%identical to) SEQ ID NO:84.

In some embodiments, the method further comprises measuring phosphatelevels in the plant.

In one embodiment, the method further comprises contacting the plantwith phosphorus oxyanions if the plant is infected with Ca. L. Thephosphorus oxyanions can be present in a solution.

In another aspect, the disclosure provides a method for treating a planthaving HLB disease. In one embodiment, the method comprises contacting aplant having HLB disease or infected with Ca. L. with phosphite orphosphorus oxyanions in sufficient amount to ameliorate the symptoms ofHLB disease in the plant. In one embodiment, the phosphite or phosphorusoxyanions are in a solution. Contacting the infected plant withphosphorus oxyanions stimulates phosphate uptake in the plant bysuppressing UBC mRNA or UBC protein expression and/or activity, which inturn induces the expression of phosphate transporters in the plant,thereby ameliorating symptoms of HLB disease.

In another aspect, the disclosure provides a kit for detection of HLBdisease or Ca. L. infection. In some embodiments, the kit comprises oneor more agents that specifically detects a small RNA, includingmiRNA399, miRNA159, siRNA1005, siRNA1008 or siRNA1009.

Definitions

As used herein and in the appended claims, the singular forms “a,”“and,” and “the” include plural referents unless the context clearlydictates otherwise. Thus, for example, reference to “a plant” includes aplurality of such plants and reference to “the tree” includes referenceto one or more trees known to those skilled in the art, and so forth.

Unless defined otherwise, all technical and scientific terms used hereinhave the same meaning as commonly understood to one of ordinary skill inthe art to which this disclosure belongs. Although any methods andreagents similar or equivalent to those described herein can be used inthe practice of the disclosed methods and compositions, the exemplarymethods and materials are now described.

Also, the use of “or” means “and/or” unless stated otherwise. Similarly,“comprise,” “comprises,” “comprising” “include,” “includes,” and“including” are interchangeable and not intended to be limiting.

It is to be further understood that where descriptions of variousembodiments use the term “comprising,” those skilled in the art wouldunderstand that in some specific instances, an embodiment can bealternatively described using language “consisting essentially of” or“consisting of.”

Two nucleic acid sequences or polypeptides are said to be “identical” ifthe sequence of nucleotides or amino acid residues, respectively, in thetwo sequences is the same when aligned for maximum correspondence asdescribed below. The terms “identical” or percent “identity,” in thecontext of two or more nucleic acids or polypeptide sequences, refer totwo or more sequences or subsequences that are the same or have aspecified percentage of amino acid residues or nucleotides that are thesame, when compared and aligned for maximum correspondence over acomparison window, as measured using one of the following sequencecomparison algorithms or by manual alignment and visual inspection. Whenpercentage of sequence identity is used in reference to proteins orpeptides, it is recognized that residue positions that are not identicaloften differ by conservative amino acid substitutions, where amino acidsresidues are substituted for other amino acid residues with similarchemical properties (e.g., charge or hydrophobicity) and therefore donot change the functional properties of the molecule. Where sequencesdiffer in conservative substitutions, the percent sequence identity maybe adjusted upwards to correct for the conservative nature of thesubstitution. Means for making this adjustment are well known to thoseof skill in the art. Typically this involves scoring a conservativesubstitution as a partial rather than a full mismatch, therebyincreasing the percentage sequence identity. Thus, for example, where anidentical amino acid is given a score of 1 and a non-conservativesubstitution is given a score of zero, a conservative substitution isgiven a score between zero and 1. The scoring of conservativesubstitutions is calculated according to, e.g., the algorithm of Meyers& Miller, Computer Applic. Biol. Sci. 4:11-17 (1988) e.g., asimplemented in the program PC/GENE (Intelligenetics, Mountain View,Calif., USA).

The phrase “substantially identical,” used in the context of two nucleicacids or polypeptides, refers to a sequence that has at least 60%sequence identity with a reference sequence. Alternatively, percentidentity can be any integer from 60% to 100%. Some embodiments includeat least: 60%, 65%, 70%, 75%, 80%, 85%, 90%, 91%, 92%, 93%, 94% 95%,96%, 97%, 98%, or 99%, compared to a reference sequence using theprograms described herein; preferably BLAST using standard parameters,as described below.

For sequence comparison, typically one sequence acts as a referencesequence, to which test sequences are compared. When using a sequencecomparison algorithm, test and reference sequences are entered into acomputer, subsequence coordinates are designated, if necessary, andsequence algorithm program parameters are designated. Default programparameters can be used, or alternative parameters can be designated. Thesequence comparison algorithm then calculates the percent sequenceidentities for the test sequences relative to the reference sequence,based on the program parameters.

A “comparison window”, as used herein, includes reference to a segmentof any one of the number of contiguous positions selected from the groupconsisting of from 20 to 600, usually about 50 to about 200, moreusually about 100 to about 150 in which a sequence may be compared to areference sequence of the same number of contiguous positions after thetwo sequences are optimally aligned. Methods of alignment of sequencesfor comparison are well-known in the art. Optimal alignment of sequencesfor comparison can be conducted, e.g., by the local homology algorithmof Smith & Waterman, Adv. Appl. Math. 2:482 (1981), by the homologyalignment algorithm of Needleman & Wunsch, J. Mol. Biol. 48:443 (1970),by the search for similarity method of Pearson & Lipman, Proc. Nat'l.Acad. Sci. USA 85:2444 (1988), by computerized implementations of thesealgorithms (GAP, BESTFIT, FASTA, and TFASTA in the Wisconsin GeneticsSoftware Package, Genetics Computer Group, 575 Science Dr., Madison,Wis.), or by manual alignment and visual inspection.

Algorithms that are suitable for determining percent sequence identityand sequence similarity are the BLAST and BLAST 2.0 algorithms, whichare described in Altschul et al. (1990) J. Mol. Biol. 215: 403-410 andAltschul et al. (1977) Nucleic Acids Res. 25: 3389-3402, respectively.Software for performing BLAST analyses is publicly available through theNational Center for Biotechnology Information (NCBI) web site. Thealgorithm involves first identifying high scoring sequence pairs (HSPs)by identifying short words of length W in the query sequence, whicheither match or satisfy some positive-valued threshold score T whenaligned with a word of the same length in a database sequence. T isreferred to as the neighborhood word score threshold (Altschul et al,supra). These initial neighborhood word hits acts as seeds forinitiating searches to find longer HSPs containing them. The word hitsare then extended in both directions along each sequence for as far asthe cumulative alignment score can be increased. Cumulative scores arecalculated using, for nucleotide sequences, the parameters M (rewardscore for a pair of matching residues; always >0) and N (penalty scorefor mismatching residues; always <0). For amino acid sequences, ascoring matrix is used to calculate the cumulative score. Extension ofthe word hits in each direction are halted when: the cumulativealignment score falls off by the quantity X from its maximum achievedvalue; the cumulative score goes to zero or below, due to theaccumulation of one or more negative-scoring residue alignments; or theend of either sequence is reached. The BLAST algorithm parameters W, T,and X determine the sensitivity and speed of the alignment. The BLASTNprogram (for nucleotide sequences) uses as defaults a word size (W) of28, an expectation (E) of 10, M=1, N=−2, and a comparison of bothstrands. For amino acid sequences, the BLASTP program uses as defaults aword size (W) of 3, an expectation (E) of 10, and the BLOSUM62 scoringmatrix (see Henikoff & Henikoff, Proc. Natl. Acad. Sci. USA 89:10915(1989)).

The BLAST algorithm also performs a statistical analysis of thesimilarity between two sequences (see, e.g., Karlin & Altschul, Proc.Nat'l. Acad. Sci. USA 90:5873-5787 (1993)). One measure of similarityprovided by the BLAST algorithm is the smallest sum probability (P(N)),which provides an indication of the probability by which a match betweentwo nucleotide or amino acid sequences would occur by chance. Forexample, a nucleic acid is considered similar to a reference sequence ifthe smallest sum probability in a comparison of the test nucleic acid tothe reference nucleic acid is less than about 0.01, more preferably lessthan about 10⁻⁵, and most preferably less than about 10⁻²⁰.

“Conservatively modified variants” applies to both amino acid andnucleic acid sequences. With respect to particular nucleic acidsequences, conservatively modified variants refers to those nucleicacids which encode identical or essentially identical amino acidsequences, or where the nucleic acid does not encode an amino acidsequence, to essentially identical sequences. Because of the degeneracyof the genetic code, a large number of functionally identical nucleicacids encode any given protein. For instance, the codons GCA, GCC, GCGand GCU all encode the amino acid alanine. Thus, at every position wherean alanine is specified by a codon, the codon can be altered to any ofthe corresponding codons described without altering the encodedpolypeptide. Such nucleic acid variations are “silent variations,” whichare one species of conservatively modified variations. Every nucleicacid sequence herein which encodes a polypeptide also describes everypossible silent variation of the nucleic acid. One of skill willrecognize that each codon in a nucleic acid (except AUG, which isordinarily the only codon for methionine, and TGG, which is the onlynatural codon for tryptophan) can be modified to yield a functionallyidentical molecule. Accordingly, each silent variation of a nucleic acidwhich encodes a polypeptide is implicit in each described sequence.

As to amino acid sequences, one of skill will recognize that individualsubstitutions, in a nucleic acid, peptide, polypeptide, or proteinsequence which alters a single amino acid or a small percentage of aminoacids in the encoded sequence is a “conservatively modified variant”where the alteration results in the substitution of an amino acid with achemically similar amino acid. Conservative substitution tablesproviding functionally similar amino acids are well known in the art.

The following six groups each contain amino acids that are conservativesubstitutions for one another:

1) Alanine (A), Serine (S), Threonine (T);

2) Aspartic acid (D), Glutamic acid (E);

3) Asparagine (N), Glutamine (Q); 4) Arginine (R), Lysine (K); 5)Isoleucine (I), Leucine (L), Methionine (M), Valine (V); and 6)Phenylalanine (F), Tyrosine (Y), Tryptophan (W).

(see, e.g., Creighton, Proteins (1984)).

The term “plant” includes whole plants, shoot vegetative organs and/orstructures (e.g., leaves, stems and tubers), roots, flowers and floralorgans (e.g., bracts, sepals, petals, stamens, carpels, anthers), ovules(including egg and central cells), seed (including zygote, embryo,endosperm, and seed coat), fruit (e.g., the mature ovary), seedlings,plant tissue (e.g., vascular tissue, ground tissue, and the like), cells(e.g., guard cells, egg cells, trichomes and the like), and progeny ofsame. The class of plants that can be used in the methods of theinvention includes angiosperms (monocotyledonous and dicotyledonousplants), gymnosperms, ferns, bryophytes, and multicellular andunicellular algae. It includes plants of a variety of ploidy levels,including aneuploid, polyploid, diploid, haploid, and hemizygous.

As used herein, the term “transgenic” describes a non-naturallyoccurring plant that contains a genome modified by man, wherein theplant includes in its genome an exogenous nucleic acid molecule, whichcan be derived from the same or a different plant species. The exogenousnucleic acid molecule can be a gene regulatory element such as apromoter, enhancer, or other regulatory element, or can contain a codingsequence, which can be linked to a heterologous gene regulatory element.Transgenic plants that arise from sexual cross or by selfing aredescendants of such a plant and are also considered “transgenic.”.

The term “promoter,” as used herein, refers to a polynucleotide sequencecapable of driving transcription of a coding sequence in a cell. Thus,promoters used in the polynucleotide constructs of the invention includecis-acting transcriptional control elements and regulatory sequencesthat are involved in regulating or modulating the timing and/or rate oftranscription of a gene. For example, a promoter can be a cis-actingtranscriptional control element, including an enhancer, a promoter, atranscription terminator, an origin of replication, a chromosomalintegration sequence, 5′ and 3′ untranslated regions, or an intronicsequence, which are involved in transcriptional regulation. Thesecis-acting sequences typically interact with proteins or otherbiomolecules to carry out (turn on/off, regulate, modulate, etc.) genetranscription. A “plant promoter” is a promoter capable of initiatingtranscription in plant cells. A “constitutive promoter” is one that iscapable of initiating transcription in nearly all tissue types, whereasa “tissue-specific promoter” initiates transcription only in one or afew particular tissue types.

A polynucleotide sequence is “heterologous” to an organism or a secondpolynucleotide sequence if it originates from a foreign species, or, iffrom the same species, is modified from its original form. For example,when a promoter is said to be operably linked to a heterologous codingsequence, it means that the coding sequence is derived from one specieswhereas the promoter sequence is derived another, different species; or,if both are derived from the same species, the coding sequence is notnaturally associated with the promoter (e.g., is a geneticallyengineered coding sequence, e.g., from a different gene in the samespecies, or an allele from a different ecotype or variety).

An “expression cassette” refers to a nucleic acid construct that, whenintroduced into a host cell, results in transcription and/or translationof an RNA or polypeptide, respectively. In the case of both expressionof transgenes and suppression of endogenous genes or regulatory elements(e.g., by small RNAs) one of skill will recognize that the insertedpolynucleotide sequence need not be identical, but may be onlysubstantially identical to a sequence of the gene from which it wasderived. As explained herein, these substantially identical variants arespecifically covered by reference to a specific nucleic acid sequence.

BRIEF DESCRIPTION OF THE DRAWINGS

FIG. 1 shows expression levels of some conserved citrus miRNAs areaffected by Las infection. Relative expression levels of some conservedcitrus miRNAs were examined at 10 wpi (a) and 14 wpi (b). Expressionlevels are presented as reads in Las-treated samples over correspondinguntreated samples. “+1” indicates a 2-fold induction; “−1” indicates a2-fold reduction; HLB: Las-infected; wpi: weeks post inoculation. Redbar highlights miR399, which is involved in citrus Pi accumulation. ThemiRNA IDs are listed at the left side.

FIG. 2 is a graph showing Ca. L. asiaticus-induced miR399 down-regulatesits target, a E2-conjugating enzyme gene UBC. UBC mRNA was measured byreal-time RT-PCR and actin was used as an internal control. Similarresults were obtained from two biological replicates.

FIG. 3 shows that Las infection causes citrus phosphorus and irondeficiency. 15-20 leaves from both untreated and Las-treated trees werecollected. For the treated plants, both asymptomatic and symptomaticleaves were collected. For symptomatic leaves a combination of blotchymottled and small chlorotic leaves of different ages was collected. Forasymptomatic and control leaves, a combination of older and youngerleaves was collected (if possible). Amount of iron, phosphorus,potassium, zinc, and carbon was examined (mean+SE). Amount of eachelement in untreated samples was assigned to 1. Experiments wererepeated for 2 times with similar results.

FIG. 4 shows Las infection down-regulates csi-PHO2 but up-regulatesPhosphate transporters (PTs). Expression levels of csi-PHO2 (UGID:1423690), csi-PT2 (UGID: 3374895), and csi-PHT2;1(UGID: 2916608)relative to citrus actin were determined by quantitative real-time PCR(mean±SE). Total RNA from both 10- and 14-wpi samples (untreated andLas-treated) was used for reverse transcription followed by quantitativePCR. Expression level of untreated samples is assigned to 1. Experimentswere repeated for 3 times with similar results. HLB: Las-infected; wpi:weeks post inoculation.

FIG. 5 shows the regulatory circuit mediated by miR399 in response toCa. L. asiaticus infection.

DETAILED DESCRIPTION OF THE INVENTION Introduction

Citrus greening or “Huanglongbing” (HLB), caused by bacteria CandidatusLiberobacter, is one of the most destructive diseases of citrus.Candidatus Liberibacter (Ca. Liberibacter or Ca. L.) is a Gram negativebacterial pathogen restricted to the phloem. The uneven distributionwithin trees and the latency of detectable symptoms make detection andconfirmation of asymptomatic infections very difficult. Therefore,developing early diagnosis biomarkers and effective reagents is anurgent need for the citrus industry, especially for those in thethreatened but un-infected regions, such as California. The recentdetection of Psyllids (the insect vector for pathogen infection) at theCalifornia-Mexico border underlines the importance of the currentinvention.

To prevent its further spread, early diagnosis before the appearance ofthe dreaded symptoms is particularly important. However, theunculturable nature of the bacteria and their low concentration anduneven distribution in the hosts make it extremely difficult to detectCa. L. infection. Instead of focusing on the bacteria, the disclosureprovides a method of early diagnosis by taking advantage of host rapiddefense responses to identify unique host biomarkers.

Some host small RNAs are rapidly and specifically induced by pathogens,which makes them one of the most attractive markers for early diagnosis.To identify HLB specific small RNAs, endogenous small RNAs were profiledby high-throughput sequencing of small RNA libraries prepared from HLBpositive and uninfected control plants. Those small RNAs that areinduced by HLB infection have the potential to serve as early detectiondiagnosis markers for HLB.

The disclosure also provides the identification of a major cause for thepathogenesis of HLB—phosphorus starvation—from the study of oneHLB-induced miRNA, microRNA399 (miRNA399), which is one of the potentialearly diagnosis markers. miRNA399 is induced specifically in response toCa. Liberibacter infection, but not to other bacterial pathogen, such asSpiroplasma citri, the causal agent for citrus stubborn disease that hassimilar symptoms. miRNA399 induction is an important response tophosphorous starvation. The level of phosphorus in Ca. L.-infectedplants was only 60-70% of that in uninfected plants. Phosphorusapplication can largely rescue the symptoms and increase the yield offruits.

Accordingly, the disclosure provides methods and compositions fortreating Ca. L. infection comprising contacting the plant with aphosphorus oxyanion solution. The disclosure also provides methods fortreating Ca. L. infection comprising downregulating expression of theubiquitin-conjugating enzyme (UBC) gene, which subsequently derepressesphosphorus transporters and increases phosphorus uptake.

By profiling small RNAs from untreated and Ca. Liberibacter-infectedCitrus sinensis (sweet orange) samples using high throughput deepsequencing, some Huanglongbing-induced small interfering RNAs (siRNAs)and microRNAs (miRNAs) were identified, providing a method of earlydiagnosis for infection. Furthermore, one of these markers, miR399, isinduced by phosphorus starvation. It targets ubiquitin E2 conjugatingenzyme genes involved in phosphorus uptake, which has been demonstratedin Arabidopsis thaliana. This result indicates that phosphorusstarvation is one of the causes for HLB disease symptom. Indeed, amarkedly decreased level of phosphorus in Ca. Liberibacter infectedcitrus was observed. Ca. Liberibacter infection causes phosphorusstarvation, and subsequently leads to quick induction of miR399, whichin turn silences UBC and derepresses the phosphorus transporters andfacilitates phosphorus accumulation. Thus, the disclosure also providesa method of treating HLB by applying phosphorus oxyanion solutions(including phosphite and phosphate salt solutions with polymeric forms)to Ca. Liberibacter-infected plants to reduce HLB symptoms.

The disclosure provides a method of determining infection of a plant ortree (e.g., a citrus plant or tree) comprising determining the level ofmiR399, siR1005, siR1008, and siR1009 in a sample from the plant ortree, wherein an increased level compared to an uninfected or healthyplant indicates that the plant or tree may be infected by Ca.Liberibacter.

The disclosure also provides a method of slowing or inhibiting thespread, growth or infection of a tree or plant with Ca. Liberibactercomprising contacting a plant with HLB disease or infected with Ca. L.with phosphorus oxyanions (which can be in solution), stimulatingphosphorus uptake in the plant by suppressing expression of UBC mRNA orUBC polypeptide in the plant, thereby ameliorating symptoms of HLBdisease in the plant.

Detection of Infection or Disease

The present methods detect the expression of small RNAs in citrus plantsthat are induced by infection with bacterial pathogens. In oneembodiment, the bacterial pathogen is Ca. Liberibacter. In oneembodiment, the bacterial pathogen is Ca. Liberibacter asiaticus (Las).

In some embodiments, infection of a citrus plant by Ca. L. induces orincreases the expression level of some small RNAs. In some embodiments,infection of a citrus plant by Ca. L. induces or increases theexpression of miRNA399, siRNA1005, siRNA1008 and/or siRNA1009 comparedto a control or healthy plant.

In some embodiments, the small RNA molecules include the sequences shownin Table 1.

TABLE 1 Sequences of small RNA molecules that areincreased following infection with Ca. L. SEQ small RNA Sequence ID NO:csi-miR399a: UGCCAAAGGAGAUUUGCCCGG  9 csi-miR399b: UGCCAAAGGAGAGUUGCCCUA10 csi-miR399c: UGCCAAAGGAGAAUUGCCCUG 11 csi-miR399d:UGCCAAAGGAGAGUUGCCCUG 12 csi-miR399e: UGCCAAAGGAGAAUUGCCCUG 11csi-siR1005: ATAGATAATGGATCAACGGTTATA 13 csi-siR1008:TCGAACAAGGTAAGGATGTCA 14 csi-siR1009: CTTCTAATAAACATGCATGAA 15

As is well known in the art, in RNA uracil (U) replaces thymide (T).Thus, T can be represented by T or U, depending if the sequence is DNAor RNA.

As is understood in the art, the sequence of siRNAs can vary at the 5′and 3′ ends. Thus, for the detection of siR1008 and siR1009, thesequence detected can be shifted either 5′ or 3′ relative to SEQ ID NOs:14 and 15. Thus, for siR1008, the detected sequence can include asequence from CCTTGTTCGAACAAGGTAAGGATGTCATTCTTT (SEQ ID NO:100), whereSEQ ID NO:14 is underlined. Likewise, for siR1009, the detected sequencecan include a sequence from CGTCTTCTAATAAACATGCATGAACTTATT (SEQ IDNO:101), where SEQ ID NO:15 is underlined.

In some embodiments, infection of a citrus plant by Ca. L. reduces ordecreases the expression level of some small RNAs. In some embodiments,infection of a citrus plant by Ca. L. reduces or decreases theexpression level of a small RNA selected from miR408, miR171, miR396,miR398, miR160, and/or miR394 compared to a control or healthy plant.

FIG. 1 shows miRNAs that have increased expression levels in citrusplants infected with Las. Thus, in some embodiments, infection of acitrus plant by Ca. L. induces or increases the expression level of amicro RNA of Table 2.

TABLE 2 miRNAs showing a relative increase in expression in citrusplants infected by Ca. L. miR399 miR159 miR393

FIG. 1 also shows that infection of a citrus plant by Ca. L. reduces ordecreases the expression level of other micro RNAs. Thus, in someembodiments, infection of a citrus plant by Ca. L. reduces or decreasesthe expression level of a microRNA shown in Table 3. In someembodiments, infection of a citrus plant by Ca. L. reduces or decreasesthe expression level of a small RNA selected from miR408, miR171,miR396, miR398, miR160, and/or miR394.

TABLE 3 miRNAs showing relative decrease in expression in citrus plantsinfected by Ca. L. miR160 miR164 miR398 miR394 miR403 miR396 miR171miR408 miR172

As is well understood in the art, miRNA molecules are thought to degradetarget mRNA or inhibit translation of RNA by precise or imprecisebase-pairing with their target RNA molecules. Thus, the sequences ofmiRNA molecules can vary slightly from their target sequence and stillfunction to inhibit gene expression. Further, the sequences of miRNA andsiRNA molecules can vary at either the 5′ or 3′ end by the addition orsubtraction of one or two nucleotides and still function to inhibitexpression of a given target.

Thus, the methods encompass detecting miRNA and siRNA molecules thathave slight variations in the nucleic acid sequence.

For example, as shown in Table 4, different sequences were identifiedfor miRNA399 in both uninfected control and Ca. L. infected citrusplants. The miR399 sequences shown in Table 4 are capable of targeting(base-pairing with) the mRNA that encodes the citrus UBC protein(accession number EY742134) or an mRNA substantially identical to themRNA that encodes accession number EY742134. In one embodiment, themiR399 sequence targets a sequence that is substantially identical toSEQ ID NO:84.

TABLE 4Citrus miRNA399 sequences and abundance upon Ca. L. asiaticus infection.SEQ Alignment Citrus smallRNA sequence ID NO: Reads subject(reverse-complemented if orientation is -) Untreated  2  2 ath-miR399bTGCCAAAGGAGAGTTGCCCTA 10 wpi  1  1 ath-miR399b TGCCAAAGGAGAGTTGCCCTG HLB 1 12 ath-miR399b TGCCAAAGGAGAGTTGCCCTG 10 wpi  2 10 ath-miR399bTGCCAAAGGAGAGTTGCCCTA  3  1 ath-miR399c TGCCAAAGGAGCGTTGCCCTG  4  1ath-miR199b TGCCAAAGGAGAGTTGCCATG  5  1 ath-miR399cTGTCAAAGGAGAGTTGCCCTG Untreated  2  2 ath-miR399b TGCCAAAGGAGAGTTGCCCTA14 wpi  8  2 ath-miR399f TGCCAAAGGAGATTTGCCCGG  6  1 ath-miR399aTGCCAAAGGAGAATTGCCCTG HLB  2 23 ath-miR399b TGCCAAAGGAGAGTTGCCCTA 14 wpi 1  8 ath-miR399b TGCCAAAGGAGAGTTGCCCTG  7  1 ath-miR399cTGCTAAAGGAGAGTTGCCCTA 99 consensusTG C/T C/T AAAGGAG A/C G/T/A TTGCC C/A T/G A/G miR399

Thus, in some embodiments, the method comprises detecting the expressionlevel of a sequence that is substantially identical to one of the miR399sequences in Table 4. For example, in some embodiments, a sequence thatis at least 60%, 70%, 80%, 85%, 90%, or 95% identical to a sequence inTable 4 is detected. In one embodiment, the method comprises detecting asmall RNA having the miR399 consensus sequence TG X₁ X₂ AAAGGAG X₃ X₄TTGCC X₅ X₆ X₇ (SEQ ID NO:99), where X₁ is C or T, X₂ is C or T, X₃ is Aor C, X₄ is G, T, or A, X₅ is C or A, X₆ is T or G, and X₇ is A or G.

In some embodiments, infection of a citrus plant by Ca. L. can bedetected by detecting the expression level of a sequence that issubstantially identical to an siRNA listed in Table 5. For example, inone embodiment, the expression level of a sequence that is at least 80%,85%, 90%, or 95% identical to an siRNA listed in Table 5 is detected.

TABLE 5 HLB induced citrus siRNAs. reads with highest copy numberuntreated HLB H/U untreated HLB H/U (SEQ ID NO:) 10 wpi 10 wpi 10 wpi14 wpi 14 wpi 14 wpi Total TTCCAGATAGAAGGCCACTCA (42) 1.0 736.0 736.01.0 106.0 106.0   844 TTCCACCAATCGATCAGGATA (43) 1.0 233.0 233.0 1.089.0 89.0   324 GCGTATGAGGAGCCATGCATA (44) 1.0 155.0 155.0 1.0 115.0115.0   272 CTTGGATTTATGAAAGACGAA (45) 1.0 127.0 127.0 1.0 1.0 1.0   130GGCAGGGCTAGTGACTGGAGTGA (46) 1.0 117.0 117.0 1.0 37.0 37.0   156ACAGGCCGCAAACATTTTCCT (47) 1.0 87.0 87.0 42.0 143.0 3.4   273ACAGACCGCACACCTTTTCTT (48) 1.0 59.0 59.0 18.0 62.0 3.4   140ATTAGGAGCTAAAATTGTTGT (49) 1.0 45.0 45.0 24.0 111.0 4.6   181ACGAAATGTGAGTAGAGTGGACAG (50) 1.0 45.0 45.0 113.0 135.0 1.2   294TTCCAAAGGGATCGCATTGA (51) 1.0 42.0 42.0 1.0 77.0 77.0   121TCATTTAAGGGTTTCGTGTTC (52) 1.0 34.0 34.0 1.0 115.0 115.0   151ACGCTCGGACGAAGCACATAGATG (53) 1.0 29.0 29.0 58.0 40.0 0.7   128TATGGGATTTACCTCGGCAAA (54) 2.0 55.0 27.5 1.0 48.0 48.0   106TGTGTGGATGAATAAGATTTC (55) 8.0 210.0 26.3 7.0 193.0 27.6   418TATCTGGATAAAAGGCTACCC (56) 212.0 5229.0 24.7 211.0 3480.0 16.5  9132TCATGGATAAGGTCATGCATT (57) 6.0 137.0 22.8 6.0 188.0 31.3   337AAAAACTTGGAAGCGTTGGAT (58) 8.0 154.0 19.3 10.0 62.0 6.2   234TCCTGCCGGGTTGCATAATCA (59) 6.0 85.0 14.2 1.0 19.0 19.0   111ATAGATAATGGATCAACGGTTATA (60) 26.0 359.0 13.8 37.0 321.0 8.7   743TCATGGATAAGGTCATGCATC (61) 25.0 308.0 12.3 24.0 403.0 16.8   760CTGAAAGCTGAGGTTGTCCTT (62) 8.0 44.0 5.5 25.0 30.0 1.2   107AGTGTCAAAAAGAGCAATGGCGTC (63) 11.0 44.0 4.0 25.0 22.0 0.9   102AATCCTTGGATTAGGAGTGTGGAG (64) 1.0 4.0 4.0 1.0 112.0 112.0   118ATCAATAAATCAGGATTGGCGGAA (65) 82.0 288.0 3.5 71.0 216.0 3.0   657CGTTAGGGAGTCCGGAGACGT (66) 13.0 44.0 3.4 20.0 37.0 1.9   114GAATAAGACATGGAGTTGGAA (67) 18.0 51.0 2.8 48.0 64.0 1.3   181AGGAAATGGACGATACGGACGCAT (68) 70.0 176.0 2.5 91.0 1.0 0.0   338TCAAGTGAGGTTCGGTCTTTGAA (69) 29.0 69.0 2.4 22.0 26.0 1.2   146TAATCGTGGGAGACGAAGCTG (70) 2184.0 5128.0 2.3 2352.0 2587.0 1.1 12251CGAAGGTCCGAGGTCGAGGTT (71) 68.0 154.0 2.3 82.0 1.0 0.0   305AGGTTTGGGCTTGTTGCAAGTAGA (72) 27.0 61.0 2.3 45.0 25.0 0.6   158TCCGGGCGGAAGACATTGTCA (73) 49.0 109.0 2.2 37.0 76.0 2.1   271AACGGAAAGAACACAACACGG (74) 736.0 1592.0 2.2 924.0 673.0 0.7  3925TGTTAGCTTTCTCGGACGCAG (75) 24.0 48.0 2.0 16.0 37.0 2.3   125TCAAGTGAGGTTCTGTCTTTGA (76) 12.0 24.0 2.0 8.0 74.0 9.3   118AAAGCAACGATTGTATGGCCA (77) 44.0 87.0 2.0 33.0 25.0 0.8   189TCAAGTGAGGTTCGGTCTTGA (78) 146.0 260.0 1.8 116.0 253.0 2.2   775GAATGTGGAATTAAGCGCACCAAA (79) 259.0 455.0 1.8 23.0 71.0 3.1   808TCGAACAAGGTAAGGATGTCA (80) 197.0 301.0 1.5 180.0 577.0 3.2  1255CTGGATGCAACTGTGGTACGG (81) 67.0 76.0 1.1 100.0 236.0 2.4   479GGTGCTTCCGGATCTCAGGAT (82) 1.0 1.0 1.0 1.0 312.0 312.0   315CACATGGGTTAGTCGATC (83) 1.0 1.0 1.0 1.0 750.0 750.0   753 H/U = ratio ofreads from HLB diseased plants and uninfected plants.

In some embodiments, infection of a citrus plant by Ca. L. is determinedby detecting the expression of one, two, three, four, five or more ofthe small RNAs from Tables 1, 2, 3, 4 or 5. In some embodiments,infection of a citrus plant by Ca. L. is determined by detecting theexpression of one, two, three, four, five or more small RNAs that aresubstantially identical to a small RNA from Tables 1, 2, 3, 4 or 5.

Methods of detecting small RNAs are well known in the art and suchmethods can be adapted to detect miRNA399, siRNA1005, siRNA1008 and/orsiRNA1009. In some embodiments, detection can include methods involvinghybridization of nucleic acids to small RNAs by base-pairing. Examplesinclude Northern analysis, polymerase chain reaction (PCR),reverse-transcription-PCR (RT-PCR) and microarrays. In some embodiments,the small RNA is detected by DNA oligonucleotides in Northern blotanalysis. In some embodiments, the small RNA is detected using a lockednucleic acid (LNA) probe in Northern blot analysis.

In some embodiments, the amount of a small RNA detected is compared tothe amount of small RNA detected in a control, uninfected plant todetermine if the plant is infected by Ca. L. In some embodiments, theamount of a small RNA detected is compared to a reference value thatcorresponds to or is indicative of the level of expression of theparticular small RNA in a plant that is not infected by Ca. L. In someembodiments, the expression of one, two, three, four, five or more smallRNAs is detected.

In some embodiments, the small RNA is detected using a nucleic acid“dipstick” or other rapid detection device.

In some embodiments, the small RNA is detected by sequencing theisolated RNA. Thus, in some embodiments, the small RNA is cloned andsequenced to determine the nucleic acid sequence of the small RNA,thereby detecting the expression of the small RNA in a plant. In someembodiments, the expression of a small RNA is detected by determiningthe number of sequencing reads that correspond to the individual smallRNA sequence. In one embodiment, the number of sequencing reads thatcorrespond to an individual small RNA sequence is compared in infectedand uninfected control plants, thereby providing an indication that theinfection increased or decreased expression of the small RNA.

In some embodiments, the small RNAs were cloned and sequenced asdescribed in the Examples. For example, total RNA is isolated from aplant tissue, 18-28 nucleotide fragments are recovered, and the purifiedfragments are ligated to adaptor oligonucleotides at the 5′ and 3′ ends.In one embodiment, the adaptor oligonucleotides serve as binding sitesfor PCR primers. The RNA fragments with the adaptor oligonucleotides arereverse transcribed and amplified by PCR. The PCR amplified products aresequenced to detect the small RNA molecules expressed by the plant.

The sequencing results revealed that the small RNAs induced by Lasinfection were from about 18 to about 28 nucleotides in length. In someembodiments, the small RNAs were identified as microRNAs (miRNAs) byaligning the sequence with conserved miRNAs in other plant species (lessthan or equal to 2 mismatches to a conserved miRNA). In someembodiments, the small RNAs were identified as miRNAs based on whethertheir precursor RNAs can form stem-loop structures. In some embodiments,the small RNA was identified as a putative small interfering RNA (siRNA)because it did not match any conserved plant miRNA or its precursor RNAdid not form stem-loop structures.

In some embodiments, the small RNA comprises a 21, 22, 23, or 24nucleotide species.

In some embodiments, increased expression is detected when theexpression level of a small RNA is at least 10%, 20%, 50%, 100%, 500% or1000% or more than the expression level detected in a control oruninfected plant. In some embodiments, decreased expression is detectedwhen the expression level of a small RNA is less than 10%, 20%, 50%,100%, 500% or 1000% or less than the expression level detected in acontrol or uninfected plant. In some embodiments, the control levelcorresponds to the expression level of a particular small RNA in ahealthy plant that does not have HLB disease. In some embodiments, thecontrol level is a reference value or average value that corresponds toor is indicative of the level of expression of the particular small RNAin a plant or population of plants that is not infected by Ca. L.

Samples

In some embodiments, the expression of a small RNA is detected in abiological sample from a citrus plant. For example, the biologicalsample can comprise bark or a leaf from an infected or control plant.

Methods of Treating

The present disclosure also provides methods of ameliorating thesymptoms of HLB infection of citrus plants. The methods described hereincan be used to reduce symptoms caused by HLB infection, includingyellowing of leaves, blotchy mottle of the leaves, zinc-deficiency-likemottle, severe chlorosis, and reduced fruit yield. It will be understoodthat symptoms of HLB vary according to the time of infection, stage ofthe disease, tree species, and tree maturity, among other things. Itwill be further understood that the disclosed methods do not necessarilyresult in eradication or cure of the infection, but can significantlyreduce the symptoms caused by HLB infection.

Thus, in some embodiments, the methods provided herein reduce thesymptoms of HLB by reducing the yellowing of leaves, resulting in agreener appearance, increasing the growth rate of the plant, and/orincreasing the fruit yield of the plant. Thus, in some embodiments, thefruit yield is improved by 5%, 10%, 20%, 30%, 40%, 50%, 60% 70%, 80%,90%, 100%, 200%, 500% of more compared to a plant that is not treatedaccording to the methods. In some embodiments, the fruit yield isincreased to 10%, 20%, 30%, 40%, 50%, 60%, 70%, 80%, 90%, or 100% of theyield of a similar plant that was not infected by HLB.

In some embodiments, the methods involve contacting an infected plantwith inorganic phosphite and/or phosphate (Pi). In one embodiment, theinorganic phosphite and/or phosphate is in a solution. In someembodiments, the plant is contacted with a phosphorus solutioncomprising phosphorus oxyanion solutions (including phosphate andphosphate salt solutions). In one embodiment, the phosphorus solutioncomprises mixed mono- and dipotassium salts of phosphorus acid, withpolymeric forms. In one embodiment, the inorganic phosphorus solution isapplied by foliar spray.

In some embodiments, the method of treatment can include transforming aplant with an expression cassette that expresses a small RNA describedherein. For example, in one embodiment, a citrus plant can betransformed with an nucleic acid that expresses miR399. While not beingbound by theory, it is believed that expression of miR399 regulates aconserved phosphate homoestasis pathway, which results in increasedtransport of Pi into the plant cells. Thus, in one embodiment, themethod comprises overexpressing miR399 in a plant. In one embodiment,the method comprises stimulating phosphate fixation or accumulation inthe plant.

In some embodiments, the method of treatment includes suppressing ordecreasing the expression or activity of a ubiquitin conjugating enzyme(UBC), for example ubiquitin-conjugating enzyme E2. In one embodiment,the citrus ubiquitin conjugating enzyme has accession number EY740382.In one embodiment, the ubiquitin conjugating enzyme is a citrus homologof the Arabidopsis PHO2 (phosphate 2) ubiquitin conjugating enzyme E2.In one embodiment, the citrus UBC is encoded by the sequence shown inSEQ ID NO:84.

In one embodiment, the method of treatment includes increasing theexpression or activity of Pi transporters. For example, in someembodiments, the Pi transporters are citrus homologs of the Arabidopsisthaliana phosphate transporter 2 (AtPT2) and phosphate transporter 2;1(AtPHT2:1). In some embodiments, the citrus phosphate transporterhomologs are represented by Unigene Csi: 14938 (UGID: 3374895) andUnigene Csi: 9842 (UGID: 2916608).

Expression Cassettes

In some embodiments, the present invention provides for expressioncassettes comprising a promoter operably linked to a polynucleotideencoding a small RNA of the invention (e.g., as described herein),wherein introduction of the expression cassette into a plant results inthe plant expressing a small RNA as described herein. In someembodiments, the promoter is heterologous to the polynucleotide. In someembodiments, the promoter is inducible. In some embodiments, thepromoter is tissue-specific.

In some embodiments, introduction of the expression cassette into aplant results in the plant having decreased expression or activity of aUBC as compared to a plant lacking the expression cassette. In oneembodiment, the introduction of the expression cassette into a plantresults in the plant having decreased expression or activity of thecitrus homolog of ubiquitin conjugating enzyme E2. In some embodiments,introduction of the expression cassette into a plant results in theplant having increased expression or activity of Pi transporters ascompared to a plant lacking the expression cassette. In one embodiment,introduction of the expression cassette into a plant results in theplant having increased expression or activity of the citrus homolog ofphosphate transporter 2 (AtPT2) and phosphate transporter 2;1(AtPHT2:1).

In another embodiment, the present invention provides for expressionvectors comprising an expression cassette of the invention (e.g., asdescribed herein).

Plants

In some embodiments, the plant is a citrus plant. In some embodiments,the citrus plant is an orange tree, a lemon tree, a lime tree, or agrapefruit tree. In one embodiment, the citrus plant is a navel orange,Valencia orange, sweet orange, mandarin orange, or sour orange. In oneembodiment, the citrus plant is a lemon tree. In one embodiment, thecitrus plant is a lime tree. In some embodiments, the plant is arelative of a citrus plant, such as orange jasmine, limeberry, andtrifoliate orange.

In some embodiments, the present invention provides for plants (or aplant cell, seed, flower, leaf, fruit, or other plant part from suchplants or processed food or food ingredient from such plants) comprisingan expression cassette comprising a promoter operably linked to apolynucleotide encoding a small RNA of the invention (e.g., as describedherein). In some embodiments, the plant has decreased UBC expression oractivity and/or increased expression or activity of Pi transporters.

Kits

In some embodiments, the disclosure provides kits that are useful fordetecting the expression of small RNAs in plants. For example, the kitcan include reagents that detect the presence of one or more small RNAsin a sample from a plant. In some embodiments, each reagent detects adifferent small RNA. In some embodiments, the reagent comprises a ligandthat is capable of specifically binding to or hybridizing with a smallRNA described herein. In some embodiments, the reagent is a nucleic acidthat is labeled with a probe or other moiety that enables detection ofthe reagent. In some embodiments, the reagents include the nucleotidesequences in Tables 1-5 above. In some embodiments, the reagents arecapable of detecting one, two, three, four, five or more of the smallRNAs described herein.

In some embodiments, the kits include primers that are useful foramplifying the small RNAs detected by the kits. For example, the kitscan include the following primer sequences:

miR399b (SEQ ID NO: 2): TGCCAAAGGAGAGTTGCCCTAmiR399b RT primer (SEQ ID NO: 85):GTCGTATCCAGTGCAGGGTCCGAGGTATTCGCACTGGATACGACATGGGC (50 nt) PCR primer:miR399b RT-F (SEQ ID NO: 86): GCGGCGGTGCCAAAGGAGAGTTmiR399b RT-R (SEQ ID NO: 87): GTGCAGGGTCCGAGGTcsi-siR1005 (SEQ ID NO: 13): ATAGATAATGGATCAACGGTTATAcsi-siR1005 RT primer (SEQ ID NO: 88):GTCGTATCCAGTGCAGGGTCCGAGGTATTCGCACTGGATACGACGTTATA (50 nt) PCR primer:csi-siR1005 RT-F (SEQ ID NO: 89): GCGGCGGATAGATAATGGATCAACGcsi-siR1005 RT-R (SEQ ID NO: 90): GTGCAGGGTCCGAGGTcsi-siR1008 (SEQ ID NO: 14): TCGAACAAGGTAAGGATGTCAcsi-siR1008 RT primer (SEQ ID NO: 91):GTCGTATCCAGTGCAGGGTCCGAGGTATTCGCACTGGATACGACATGTCA (50 nt) PCR primer:csi-siR1008 RT-F (SEQ ID NO: 92): GCGGCGGTCGAACAAGGTAAGGcsi-siR1008 RT-R (SEQ ID NO: 93): GTGCAGGGTCCGAGGTcsi-siR1009 (SEQ ID NO: 15): CTTCTAATAAACATGCATGAAcsi-siR1009 RT primer (SEQ ID NO: 94):GTCGTATCCAGTGCAGGGTCCGAGGTATTCGCACTGGATACGACCATGAA (50 nt) PCR primer:csi-siR1009 RT-F (SEQ ID NO: 95): GTGCAGGGTCCGAGGT.

EXAMPLES Example 1

This example describes methods for detecting and treating citrus plantsinfected with Ca. L. asiaticus (“Las”).

Materials and Methods

Plant Material Maintains, Inoculation, and RNA Extraction

For examining Las induced small RNA, two-year-old greenhouse-grown‘Navel’ orange (C. sinensis) scions on Cleopatra mandarin (C.reticulata) rootstocks were inoculated by grafting a combination of treebark pieces or leaf pieces onto the rootstock portion of each plant.Five plants were inoculated with bacterium-free tissue pieces and 19plants were inoculated with infected tissue for 10 and 14 weeks postinoculation (wpi), respectively. For infected plants, bark and leafpieces were obtained from infected greenhouse-grown Valencia scions. Fornon-infected plants, bark and leaf pieces were obtained from a healthygreenhouse-grown Valencia orange, confirmed negative for Las by PCR. Forexamining Spiroplasma citri induced small RNA, Spring Navel/Carrizo wasgrafted to Madam Vinous (C. sinensis) receptor plants, which arepositive for S. citri by culturing and PCR using primers to spiralingene sequences. Before inoculation plants were fertilized using a17N-6P-10K controlled release fertilizer (Scotts Sierra, Marysville,Ohio, USA). Plants were arranged randomly on the greenhouse bench andkept under natural light conditions at a temperature of 17-25° C.

Small RNA Library Construction

Total RNA was resolved by a denaturing 14% polyacrylamide gel. 18-28 ntRNA fragments were recovered. The purified fragments were 5′ adapted(GUU CAG AGU UCU ACA GUC CGA CGA UCA G; SEQ ID NO:97), gel purified, and3′ adapted (UC GUA UGC CGU CUU CUG CUU G; SEQ ID NO:98). The RNAfragments adapted at both ends were gel purified and reverse transcribed(SuperScript II, Invitrogen). After PCR amplification for 15 cycles (98°C. for 30 sec; 98° C. for 10 sec; 60° C. for 30 sec), the PCR productswere gel purified and sequenced according to Solexa small RNA sequencingprotocol.

Mineral Measurement

For mineral analysis, six to eight leaves were collected from each plantfor analysis. Leaves from two non-infected plants of the same age andthe same scion/rootstock combination were collected for comparison. Leaftissue was washed for 15 sec each in RO-D (reverse osmosis-distilled)water, 0.01% detergent (Citranox, Alconox, Inc., White Plains, N.Y.) and0.1 N HCl solution, followed by three more rinses in RO-D water, driedat 80° C. for at least 24 h in a forced-air oven, dry weight recorded,and leaf tissue milled to pass a 20-mesh screen. Leaf tissue (500 mg)was digested in 10 ml of concentrated HNO3 (trace metal grade) at 300psi and 170° C. for 10 min in a microwave (model Mars 5, CEM Corp.,Mathews, N.C.). Leaf digestates were brought to volume in 100 mlvolumetric flasks and filtered (no. 41; Whatman Paper, Maidstone, Kent,U.K.). Foliar levels of phosphorous (P), potassium (K), iron (Fe),manganese (Mn), copper (Cu), zinc (Zn), sodium (Na), calcium (Ca), andmagnesium (Mg) were determined by inductively couple plasma (ICP [IRIS1000 HR Duo, ThermoElemental, Franklin, Mass.]) and nitrogen (N) andcarbon (C) by dynamic-flash combustion and GC separation (CNS analyzer[NC 2100, CE Elantech, Inc., Lakewood, N.J.]).

Application of Inorganic Phosphorus Solution

The foliar nutrient sprays started with the spring flush in 2008.Applications have been made three times each year and timed with theinitiation of the new vegetative flushes in the spring (March), summer(June), and autumn, (September). The foliar spray treatment has been a3-18-20 analysis liquid fertilizer of potassium poly phosphate with 1.0pint/gal. mono- and dipotassium salts of phosphorus acid 56% (K-Phite,Plant Food Systems, Inc., Zellwood, Fla. 32798), plus 8.4 lbs. spraygrade potassium nitrate (KNO₃), and 5.0 gallons 435 citrus spray oil in250 gal./acre rate applied by hand gun to the foliage until runoff.Untreated control trees received no foliar applied nutrients. All treesreceived two ground-applied applications of a controlled-releasefertilizer to the soil each year of a 14-0-24 analysis. The nutrientsare from sulfur coated urea, sulfur coated sulfate of potash, ammoniumnitrate, sulphate of potash magnesia, and muriate of potash. Thefertilizer contained 2.3% Mg, 005% B, and 7.57% S. Insect management wasdetermined by scouting and appropriate insecticides were applied bytractor pulled speed sprayer when required. Copper was used as afungicide for disease control. Irrigation was by micro sprayers undereach tree. Systemic and pre-emerge herbicides were used for weedmanagement. Other standard good grove management practices known tofavor good production were used.

PCR Detection of Las

The presence of Las was tested monthly and symptom-development wasrecorded. Petioles were ground in liquid nitrogen with a mortar andpestle. One hundred mg of ground tissue was used for DNA extractionaccording to the manufacturer's instructions (Qiagen, Valencia, Calif.,USA). PCR amplifications were performed as described by Albrecht andBowman (2008).

RNA Extraction

Total RNA was extracted from fully expanded leaves according to Strommeret al. (1993) with slight modifications.Phenol/chloroform/isoamylalcohol (25:24:1) extraction was followed byone extraction with chloroform/isoamylalcohol and precipitation of RNAwith isopropanol at −20° C. overnight. RNA was pelleted bycentrifugation at 10,000 g and 4° C. for 1 h. RNA pellets were stored at−80° C. and shipped to UC Riverside for construction of small RNAlibraries. Leaf petioles were used for PCR detection of the bacterium.

Small RNA Parsing

After Illumina sequencing, the generated sequences (in FASTA format)were parsed and trimmed in total. An R program was used to identify thesmall RNA/adaptor boundaries, trim off the adaptors with variant length,and sort the sequences to respective source libraries. After removingthe sequences shorter than 18nt or longer than 28nt, the remainingunique small RNAs were filtered by BLAST against a citrus non-coding RNAdatabase, which included rRNA, tRNA, snoRNA and snRNA. The citrusnon-coding RNA database was generated by combining Rfam database(Release 9.1, http://www.sanger.ac.uk/Software/Rfam/) with NCBI citrusEST sequences (http://www.ncbi.nlm.nih.gov/). The resulted small RNAsequences were interpreted as citrus small RNA and were studied in thisresearch.

Small RNA Cluster Analysis and Annotation

Citrus EST sequences from NCBI (http://www.ncbi.nlm.nih.gov/) andnon-redundant short EST sequences (<150 bps in length) parsed from thetrace files (provided by Dr. Tim Close; UC Riverside) were assembledusing CAP3 to generate a citrus EST Uniset database. The parsed smallRNA sequences were then aligned to citrus EST Uniset database withBOWTIE (version 0.10.0, http://bowtie-bio.sourceforge.net/index.shtml),allowing up to 2 mismatches. Such generated sequences were mapped, andthe reads for each source library was counted. To compare the expressionlevel between different libraries, small RNA counts were normalized tothe original library size. Cluster analysis on small RNAs matching tocitrus EST Uniset was done by R program. A small RNA cluster was definedby aligning to the same locus with at least 16nt overlapping. Theclusters featured with forward (F) and reverse/complementary (R/C)alignments were counted in the 4 libraries, respectively. The functionalannotation for citrus EST Uniset genes was from NCBI non-redundantprotein database, Arabidopsis protein database (TAIR8_pep_20080412) andrice protein database (Version 5.0). Only the Uniset genes aligned withmore than 10 small RNA copies were listed in the final small RNAexpression table.

Conserved and Novel miRNA Analysis

For identifying conserved microRNAs, candidate small RNAs were analyzedby BLAST against microRNA Registry and Plant MicroRNA database (PMRD;http://bioinformatics.cau.edu.cn/PMRD/), respectively. For predictingnovel miRNAs, each EST sequences matched with small RNAs were passedthrough an R script that detects a 450 bp pile-up matching region, andretrieve the flanking sequences by 200 bp from each side of the region.All retrieved EST fragments were folded with RNAfold (version 1.6.1).Structures of the EST sequences with minimum free energies were furtheranalyzed by using an R and Perl script to retrieve stem-loops from thesecond structures and to check whether these stem-loops satisfy thefollowing criteria: (1) the length of stem is longer than 20 bp; (2) nomore than 4 bugles in a stem; (3) no more than 3 bp mismatches in abugle; (4) putative microRNA locates on one strand of the stem, whilemiRNA* on the complementary strand; (5) no small RNA read matches to theloop region. Such predicted miRNAs were subjected to northern blotvalidation.

Conserved and Novel Citrus MicroRNA Target Prediction

MicroRNA targets were computationally predicted from the citrus ESTUniset using TargetFinder program (Release 1.5,http://jcclab.science.oregonstate.edu/node/view/56334) with defaultparameter settings. Briefly, potential targets from FASTA searches(+15/−10 match/mismatch scoring ratio, −16 gap penalty, and a RNAscoring matrix) were scored using a position-dependent, mispair penaltysystem. Penalties were assessed for mismatches, bulges, and gaps (+1 perposition) and G:U pairs (+0.5 per position). Penalties were doubled ifthe mismatch, bulge, gap, or G:U pair occurred at positions 2 to 13relative to the 5′ end of the microRNA. Only one single-nt bulge orsingle-nt gap was allowed. Based on a reference set of validatedmicroRNA targets, only predicted targets with scores of four or lesswere considered reasonable. The functional annotation for targets ofconserved and novel microRNA was done from NCBI non-redundant proteindatabase.

Northern Hybridization

Fifteen to 100 micrograms of total RNA were loaded per lane, dependingon signal strength. The RNAs were resolved on a denaturing 14%polyacrylamide gel and electro-blotted onto Hybond N⁺ membranes(Amersham) overnight in a cold room at a constant 150 milliampere(Bio-Rad). Membranes were cross-linked with 0.15M1-ethyl-3-[3-dimethylaminopropyl] carbodiimide hydrochloride (EDC) at60° C. for 2 hours and dried at 80° C. for 2 hours. Citrus miRNA probeswere synthesized as their reverse complementary sequences to the maturemiRNAs and were labeled at the 5′ end with γ⁻³²P-ATP using T4polynucleotide kinase (Biolabs). Membranes were pre-hybridized usingPerfect Hyb Plus buffer (Sigma) for 2 hours and hybridized with probesovernight at 37° C. Membranes were exposed overnight either to storagephosphor screens (GE health) or X-ray films after being washed fourtimes (20 min each).

Results

Small RNA Profiling of HLB-Positive and Healthy Sweet Orange

To study citrus small RNAs in responding to Las infection, C. sinensisplants were infected by Las by graft. Samples were collected at an earlytime point (10 wpi) and a late time point (14 wpi), respectively. SmallRNAs ranging from 18 to 28 nucleotides (nt) from both healthy andLas-infected (treated) plants were cloned and sequenced. About half ofthe sequences can be aligned to the currently available citrus genomethat was assembled from different databases. The identified citrus smallRNAs were grouped into miRNAs and siRNAs by whether the sequence can bematched to conserved miRNAs (=<2 mismatches) in other plant species (top20 plant species in miRBase; Table 6), and/or whether their precursorRNAs can form stem-loop structures. Our data showed that about 30% ofthe deep-sequencing reads belong to miRNA, while another 70% belong tosiRNA due to the lack of miRNA characteristics. Since these citrussiRNAs didn't show any sequence similarities to siRNAs identified inother plant species or other organisms, we interpret them as C. sinensisspecific siRNAs and named them as csi-siRNAs. Our deep-sequencingresults showed that, similar to the model plant such as Arabidopsis,citrus small RNAs are featured by 21 nt and 24 nt species, as well asrelatively low-abundant 22 and 23-nt populations. The majority of theidentified citrus miRNAs start with a uracil (U), while most of thesiRNAs start with both U and A (adenine).

TABLE 6 conserved citrus miRNAs and abundance in libraries. miRNA 10 wpi14 wpi families Untreated HLB Untreated HLB Total miR157 31896 3402046477 58511 170904 miR166 2260 2345 69580 58393 132578 miR167 8546 19162621 7820 20903 miR164 4409 1775 2271 1227 9682 miR172 2696 1440 23661199 7701 miR165 123 108 4620 2790 7641 miR168 1160 1827 1393 2324 6704miR159 109 749 321 1326 2505 miR396 707 113 1173 139 2132 miR845 855 234311 209 1609 miR169 220 84 343 424 1071 miR403 612 104 104 42 862 miR827223 124 189 156 692 miR160 159 66 406 41 672 miR156 86 43 236 283 648miR170 215 129 97 195 636 miR162 60 49 161 147 417 miR393 16 66 22 271375 miR171 191 22 39 15 267 miR394 23 8 101 19 151 miR399 4 26 6 33 69miR398 10 4 23 5 42 miR408 14 1 7 3 25 miR390 7 4 3 4 18 miR395 2 6 0 513 miR397 3 0 3 4 10

Citrus miRNAs were searched against database containing miRNA sequencesfrom Arabidopsis thaliana, Brassica rapa, Populus trichocarpa (westernbalsam poplar), Gossypium hirsutum (upland cotton), Brassica napus,Glycine max, Vitis vinifera (wine grape), Solanum lycopersicum, Medicagotruncatula, Oryza sativa, Zea mays (maize), Triticum aestivum (breadwheat), Sorghum bicolor, Saccharum officinarum (noble cane), Hordeumvulgare subsp. vulgare, Pinus taeda, Ricinus communis, Selaginellamoellendorffii, Physcomitrella patens, and Chlamydomonas reinhardtii.miRNAs with more than 10 reads were deemed as valid conserved miRNAs andreported in Table 6.

Research has shown that some Arabidopsis siRNAs that are specificallyelicited upon bacterial pathogen infection participate in host innateimmunity by regulating the expression level of their targets(Katiyar-Agarwal et al., 2006; Katiyar-Agarwal et al., 2007; Zhang etal., 2011). We investigated whether csi-siRNAs expression is alteredupon Las-infection. We compared csi-siRNAs expression levels in bothhealthy and Las-infected samples and found that, similar in Arabidopsis,some csi-siRNAs also showed varied expression levels by Las infection(Table 5). Csi-siRNAs identified in this research were sorted accordingto their expression levels between healthy and treated samples and 10representatives (csi-siR1001-10) were selected for experimentalvalidation. Csi-siR1005 was induced upon Las infection both at 10- and14-wpi. Csi-siR1005 can target a citrus protein without known function.Similarly, csi-siR1008 and csi-siR1009 also showed weak but noticeableinduction upon Las infection. Csi-siR1008 can target a citrus proteinhomolog to an Arabidopsis putative disease resistance protein (CC-NBSclass) with perfect match, whereas csi-siR1009 targets a putativedisease resistance protein with 2 mis-matches. It suggests thatcsi-siR1008 and csi-siR1009 may play potential roles in citrus innateimmunity by regulating the expression of disease resistance proteins.These inductions are specific to Las infection since when samples weretreated with S. citri that causes citrus stubborn disease, nosignificant induction was observed.

We focused our study on citrus miRNAs, which have been more extensivelystudied compared to siRNAs. A potential advantage of focusing on miRNAsis that currently available information in other systems may facilitatethe study in citrus. Our analysis showed that most of the citrus miRNAs(conserved plus novel miRNAs) are 21-nt long, and dominantly favored a Uat the 5′ end. Among the citrus miRNAs identified, about 240 miRNAsbelonging to at least 26 miRNA families (reads>=10) could also be foundin other plant species including Arabidopsis thaliana (Table 6),indicating that regulatory mechanisms in other plant species may alsoapply to the citrus system. Beyond these conserved miRNAs, 13 identifiedcitrus miRNAs could not be matched to any currently known plant miRNAs,and therefore were termed as novel miRNAs. Similarly, they were named ascsi-miRNAs (Table 7).

TABLE 7 Novel citrus miRNAs. Candidate microRNA sequence microRNA staruntreated HLB untreated HLB ID (SEQ ID NO:) (SEQ ID NO:) 10 wpi 10 wpi14 wpi 14 wpi Csi- UGAAGCUGCCAGCAUGA AUCAUCUGGCAGUUUC 7713.00 9946.31625.47 776.91 miR5001 UCU (16) ACC (17) Csi- UAGAUAAAGAUGAGAGAUUUUCUCUUAUCGUUA 604.00 104.13 675.33 235.24 miR5002 AAAA (18)UCUGU (19) Csi- UUUGUUGCAUGAUGCUG CUACCCGCAUCAUGCA 99.00 212.42 227.7553.65 miR5003 AUAA (20) ACAAA (21) Csi- AGUGUUAGGUGUAGAGAUCUCGUACUUCUCUUC 125.00 122.87 132.57 99.05 miR5004 AGCACGA (22)ACCAAGCA (23) Csi- AUUCGGGACGAGUUUAC CGUAAACUCGUCUCGU 83.00 97.88 100.8550.56 miR5005 AAC (24) ACUU (25) Csi- AAUGGCUGGAUCCAGCU ACAAUUGGAUUUAGCC51.00 42.69 48.72 25.79 miR5006 GUGG (26) AUUAA (27) Csi-AUGCCGUAUCACGUGGG CCCACUUGAUGUGUCA 44.00 49.98 35.13 27.86 miR5007AG (28) UUC (29) Csi- UUUGAUGCCUUCUUUAG GAUUAAGGAGAGUUUU 5.00 10.4118.13 10.32 miR5008 UCGC (30) CAGU (31) Csi- GACAGAAGAGAGTGAGCGCUCGCUCCUCUUCUG 0.00 5.21 0.00 31.98 miR5009 AC (32) UCAG (33) Csi-GUAUAUAUCUUGCAUGC GUAUGUAAGAUACAUC 12.00 9.37 5.67 2.06 miR5010 AUG (34)CCC (35) Csi- UUUCUCUUAUCGUUAUC UAGAUAAAGAUGAGAG 17.00 1.04 3.40 3.10miR5011 UGU (36) AAAAA (37) Csi- UUGUUGUUGAGUGUGUA ACAUAAAUACUUAAUA 8.008.33 2.27 2.06 miR5012 UGUUA (38) AUAAUC (39) Csi- UCGUCCUUCUCUCAUAUAGAGUGGGUGGGUGGA 1.00 4.17 6.80 1.03 miR5013 UUUU (40) GAGG (41) *Readsare normalized to library sizes. Csi = citrus sinensis. wpi = weeks postinfection.

The computationally predicted novel csi-miRNAs were experimentallyvalidated using radioactively labeled probes. The Northern hybridizationresults showed a faithful agreement with the deep-sequencing results.For example, csi-miRNA5001, 5002, 5003, 5004, 5005, 5006, and 5007,which showed moderate or high reads in our deep-sequencing analysis,could be easily detected by their corresponding probes in Northernblots. In contrast, other csi-miRNAs that were predicted as low abundantappeared to be too weak to be detected (data not shown). This wasconfirmed by increasing the detecting limit by employing a lockednucleic acid (LNA) probe: when probed by a LNA probe, csi-miRNA5009 wasclearly detected. However, there is no noticeable variation betweenhealthy and treated samples, or between 10 wpi and 14 wpi time points,indicating that under the examined conditions, these csi-miRNAs may notbe subject to expression alteration.

Some Citrus miRNAs and siRNAs are Differentially Expressed in Healthyand HLB-Positive Plants

Computational analysis of our data revealed that some of the conservedcitrus miRNAs showed elevated expression upon Las infection, whileothers showed reduced expression. For example, at 10 and 14 wpi, miRNAssuch as miR159, 399, 393 are noticeably induced (>2 fold), while miRNAssuch as miR160, 396, 394, 398, 171, 403, and 408 are clearly reduced (>2fold), as shown in FIGS. 1a and 1b . This implies that these csi-miRNAsmay be subject to the influence of Las infection, or even play somepotential roles in Las-elicited plant immunity. We tested some of themiRNAs by using Northern hybridization. Our results showed that allthese alterations in miRNA expression upon Las infection weresuccessfully validated using sequence-specific probes against the familymember with the highest reads (data not shown). This induction isspecifically due to Las infection since S. citri infection didn't shownoticeable alteration on miR399 and miR159 expression level.

Induction of miR399 Revealed Phosphorus Deficiency in HLB-PositivePlants

In spite of its relatively low abundance, miR399 showed distinguishableinduction upon Las challenge (FIG. 1). MiR399 has been shown to beinvolved in phosphate (Pi) homeostasis in Arabidopsis (Fujii et al.,2005; Lin et al., 2008; Hsieh et al., 2009). Some studies have shownthat miR399 is the phloem-mobile long-distance signal involved inphosphate starvation response (Buhtz et al., 2008; Pant et al., 2008;Buhtz et al., 2010). In Arabidopsis thaliana, miR399 has multiple targetsites in the 5′UTR of the transcript of a gene encoding a putativeubiquitin-conjugating enzyme (PHO2; At2g33770), which in turn negativelyregulates Pi transporters (PT) (Fujii et al., 2005; Bari et al., 2006;Lin et al., 2008). Upon Pi deficiency, the miR399-mediated Pihomeostasis mechanism was turned on: the Pi-deficiency-induced miR399down-regulates PHO2, which releases its inhibitory role on PTs; theincreased PTs transport more Pi into the cells as a consequence, whichalleviates Pi deficiency (Fujii et al., 2005; Lin et al., 2008).Research showed that the same Pi regulatory mechanism might also existin other plant species such as rapeseed (Buhtz et al., 2008) and pumpkin(Pant et al., 2008). Interestingly, Bari et al. (2006) identifiedpotential orthologs of PHO2 from orange (Citrus aurantium), implyingthat the regulatory mechanism may be also conserved in citrus plants.

The variation on citrus miR399 level between healthy and Las-treatedcitrus samples, and the potential conservation of the Pi regulatorymechanism prompted us to investigate the Pi level in these samples. Wehypothesize that if the same Pi-regulatory mechanism also applies in C.sinensis, then the Las-infected C. sinensis would show deficiency in Pilevel and elevated miR399 expression level; and we would also observereduced PHO2 (citrus homolog) mRNA and increase PTs.

We collected leaves from both healthy and Las-infected plants(PCR-confirmed). Our measurement showed that the phosphorus level inleaves of treated plants was more than 35% lower than in leaves fromuntreated plants. We also observed iron (Fe) deficiency in Las-infectedplants, which indicates the phloem localized bacterial pathogen conferssimilar negative influence on several mineral elements. However,homeostasis of entire mineral elements is not negatively affected sincesome elements showed no noticeable variation, such as potassium (K) andzinc (Zn) (FIG. 4a ). In contrast, copper (Cu) accumulated more inLas-infected plants (data not shown).

We further tested our hypothesis by detecting PHO2 and PTs transcriptsin both healthy and Las-infected plants by real-time PCR. After databasesearching, we identified one citrus PHO2 and two citrus PTs: csi.2677(UGID: 1423690; homolog to Arabidopsis PHO2 [At2g33770;identity=69.9%]), csi.14938 (UGID: 3374895; homolog to AtPT2[Arabidopsis thaliana phosphate transporter 2; AT2G38940;identity=85.2%]), and csi.9842 (UGID: 2916608; homolog to PHT2;1[phosphate transporter 2;1; AT3G26570; identity=88.3%]). These citrushomologs therefore are referred as csi-PHO2, csi-PT2, and csi-PHT2;1hereafter. We examined gene expression levels of PHO2 and PTs byreal-time PCR. As shown in FIG. 4, we observed reduced csi-PHO2expression level at both 10- and 14-wpi, agreeing with the increasedmiR399 expression at both time points, as well as the reduced Pi level;we also observed elevated expression of csi-PT2 and csi-PHT2;1, whichare opposite to csi-PHO2 and in agreement with csi-PHO2's repressiverole on csi-PTs. Taken together, our hypothesis was validated by theobservation in concord with all the predictions: Pi-deficiency andinduced miR399 in Las-infected plants, opposite expression profiles ofcsi-PHO2 (reduced) and csi-PT2 and -PHT2;1 (induced) at both time pointsof the Las-infected plants. Therefore, we propose that the miR399-PHO2regulatory machinery is a conserved Pi homeostasis regulatory mechanism,at least between Arabidopsis and citrus.

Application of Inorganic Phosphorus Solution Largely Reduces HLBSymptoms

If Las infection leads to host Pi deficiency, which consequently causedobserved symptoms such as yellowing, blotchy mottle,zinc-deficiency-like mottle, severe chlorosis, and most importantly,reduced fruit yield, then applying Pi to the infected plants shouldalleviate at least some of the HLB symptoms. We foliarly appliedpotassium poly phosphate (mono- and di-potassium salts of phosphorusacid 56%, plus potassium nitrate [KNO₃], and citrus spray oil) toLas-infected plants for two years. As controls, we also applied KNO₃ andcitrus spray oil to infected plants, as well as untreated control treesthat received no foliarly applied nutrients. After applying potassiumpoly phosphate for two years, the Las-infected plants showeddramatically reduced disease symptoms. Compared to untreated controltrees and trees treated with nutrients only, the phosphate-treated treeshave greener appearance, vigorous growth, and increased fruit yield. Onthe leaves from trees applied potassium poly phosphate, there is hardlyany yellowing or blotchy mottle visible. In contrast, these symptoms arequite obvious on leaves from untreated trees or trees treated withnutrients only. The observed reduced symptoms are not due to applicationof potassium since the control plants also received KNO₃. Therefore, weconclude that applying Pi, but not potassium nutrients, to Las-infectedcitrus can relieve HLB symptoms. Our results indicate that Piapplication may be employed for HLB management and restore fruit yieldin infected regions.

DISCUSSION

Citrus is an important economic plant globally, which is currentlythreatened by the very destructive citrus disease-HLB (Bove, 2006). Oneof the current challenges for HLB management is an effective HLB diseasecontrol reagent that can replace the costive eradiation procedure, whichis the only option currently available. Plants have evolved multiplelevels of immune responses, including basal defense triggered byvirulent pathogens in susceptible hosts and resistance (R) gene-mediatedresistance activated by avirulent pathogens in resistant hosts (Chisholmet al., 2006; Jones & Dangl, 2006). Infection of bacterium Ca. L. indifferent genotypes of citrus plants causes different degrees of diseaseand symptoms. Although there is no known complete resistance in Citrusspp., tolerant citrus cultivars that have a very mild or no obviousdisease symptom and with low bacterial titer were identified (Folimonovaet al., 2009). Thus citrus hosts are capable of recognizing Ca. L.infection and responding to the pathogen in different degrees. Thesignificance of this research is investigating the involvement of thecitrus innate immune system in Ca. L. infection. Studying andunderstanding how citrus hosts use innate immunity as tools against Ca.L. should lead to utilizing of these innate immune tools as efficientbiological reagents against Ca. L. in future. This is going to have ahuge impact on the citrus industry considering HLB management-relatedcost is estimated about 40% higher than pre-HLB costs in the Unitedstates (Irey et al., 2008) and about 12.65-38.73% of the totaloperational costs in Sao Paulo, Brazil (Belasque et al., 2010).

Previous studies have shown that small RNAs are involved in basaldefense and R gene-mediated resistance (Katiyar-Agarwal et al., 2006;Navarro et al., 2006; Katiyar-Agarwal et al., 2007; Zhang et al., 2011),as well as plant fitness adjusting mineral homeostasis (Jones-Rhoades &Bartel, 2004; Fujii et al., 2005; Burkhead et al., 2009). Ourinvestigation of citrus small RNA populations revealed that citruspossess many conserved miRNAs, which can also be found in other plants,such as Arabidopsis. In most of the cases, bioinformatics analysis dataagreed with our experimental very well. For example, there is a verygood match between the bioinformatics prediction and experimentalvalidation of the expression level of citrus conserved miRNAs; most ofthe miRNAs that were predicted more than 2-fold increase or decrease inboth 10- and 14-wpi (FIG. 1) were successfully experimentally validated.This indicates that our bioinformatics approaches are not only powerful,but also reliable in investigating the dynamic small RNA populationsupon pathogen challenge. Our results showed that some citrus miRNAs,such as miR 160, 396, 398, and 399, demonstrated distinguishedexpression patterns between healthy and Las-infected plants, among whichmiR399 and miR396 are the ones with potential practical value due totheir relative high abundance.

Furthermore, the Las elicited miR399 induction and Pi deficiencysuggests there might be connections between HLB and Pi level. Researchhas shown that in Arabidopsis the miR399/PHO2 regulatory machinery isturned on upon Pi starvation, and over-expression of miR399 inducesaccumulation of Pi (Fujii et al., 2005; Lin et al., 2008; Buhtz et al.,2010). Using public genomic DNA and expressed sequence tag data, Bari etal (2006) assembled potential orthologs of PHO2 from rice, Medicagotruncatula, poplar, wheat (Triticum aestivum), soybean (Glycine max),cotton (Gossypium hirsutum), apple (Malus domestica), and orange (Citrusaurantium), but not in the Physcomitrella patens EST database or thegenome sequences of P. patens and C. reinhardtii. These results suggestthat the regulatory mechanism may be conserved across angiosperms andthat it may have emerged during the evolution of higher plants.

If the citrus and Arabidopsis miR399 are really functionally conserved,as suggested by our data and observations by Bari et al (2006), this canbe interpreted as that Las infection reduces host Pi level, which inturn triggered miR399 induction, as observed in this study. If this istrue, then citrus miR399 should be induced upon Las infection, shouldparticipate in host Pi up-regulation pathways (by regulating PHO2directly, and csi-PT2 and -PHT2;1 indirectly), and applying Pi shouldalleviate HLB symptoms.

Several lines of evidence support our hypothesis. First, as wedemonstrated in this study, when citrus was affected by Las, Pideficiency was observed (FIG. 4a ), which is about >30% lower than inthe healthy leaves. Some mineral nutrients, such as zinc and potassium,did not show noticeable variations. Second, citrus miR399 wasspecifically induced upon Las (FIG. 1), but not S. citri challenge.Although the overall expression level is relatively low in citrus,elevation in expression level could be detected unambiguously andconsistently. Third, we proved that citrus miR399 is functionallyconserved with its Arabidopsis homolog, which regulates Pi accumulationthrough down-regulating csi-PHO2 (FIG. 4). We also observed elevatedexpression of two Pi transporters-csi-PT2 and -PHT2;1, respectively. Thecognate induction of citrus miR399, reduction on csi-PHO2 level, andinduction of csi-PT2 and -PHT2;1 after Las infection suggests there is afunctional miR399/PHO2/PTs pathway controlling Pi homeostasis, as beingrevealed in other species such as Arabidopsis, rapeseed, and pumpkin(Fujii et al., 2005; Buhtz et al., 2008; Pant et al., 2008). Fourth,field practice showed that when Pi is applied, HLB symptom isalleviated. This indicates that Las infection caused Pi deficiencycontributes to observed HLB symptoms, including yellowing and reducedyields, and when normal Pi level is restored by applying Pi (or byendogenous mechanisms, such as the miR399 pathway), Las-infected citruswould show reduced symptom.

Based on our results, we propose a model in which citrus miR399 may playa role in citrus response to Las infection. In this model, Las infectioncauses a Pi deficiency on the host. This is consistent with currentobservation that Las is restricted to phloem, and Las infection usuallycause phloem congestion (Folimonova & Achor, 2010). Reduced Pi levelinduces miR399 expression, which in turn down-regulates its target,csi-PHO2, a citrus ubiquitin conjugating enzyme. Based on its homologyto Arabidopsis PHO2, csi-PHO2 may be a suppressor that restricts citrusphosphate transporter activity at normal conditions. Upon Las infection,the miR399-mediated reduction of csi-PHO2 consequently leads toinduction of phosphate transporters, which recharge host cells' Pireservoir. This model is similar to the model proposed in Arabidopsis,where mR399 and PHO2 play important roles in response to Pi deficiency(Lin et al., 2008). Foliar applying Pi to the Las-infected plants, whichreduced HLB symptoms and enhanced fruit yield, validated our model. Ourmodel may explain the filed practice that applying nutrition to infectedtree can maintain tree health and productivity. Studies showed that Ca.L. infection restricted either nutrient uptake or transport and thatfoliar applied minerals could prolong tree life and reduce yield losses(Pustika et al., 2008).

Our model does not exclude the possibility that miR399 can be directlyinduced by Las infection. Research has shown that miR393 and miR393b*can be induced by bacterial pathogens, and their induction contributesto antibacterial resistance (Navarro et al., 2006; Zhang et al., 2011).We speculate that over-expression of miR399 in Las-infected citrusshould reduce HLB symptoms. Whether normal or enhanced Pi level incitrus plants have an inhibitory effect on infection and/or propagationof Las in hosts remains to be investigated in future. It should be notedthat besides phosphorus, we also observed reduced iron level after Lasinfection (FIG. 4a ). It seems that after Las infection, the congestedphloem also affects homeostasis of other mineral elements.

This example demonstrates the Ca. L. asiaticus infection of citrusplants can be detected by detecting the expression of miRNAs and siRNAs,and that HLB disease symptoms can be ameliorated by treatment of plantswith phosphorous.

It is understood that the examples and embodiments described herein arefor illustrative purposes only and that various modifications or changesin light thereof will be suggested to persons skilled in the art and areto be included within the spirit and purview of this application andscope of the appended claims. All publications, sequence accessionnumbers, patents, and patent applications cited herein are herebyincorporated by reference in their entirety for all purposes.

What is claimed is:
 1. A method for detection of Huanglongbing (HLB)disease or Ca. L. asiaticus-infection in a citrus plant, comprisingdetecting in a sample from the citrus plant the level of expression ofone or more RNA selected from the group consisting of miRNA399,siRNA1005, siRNA1008 and siRNA1009, wherein increased expression of theone or more RNA compared to expression of the RNA in a non-infectedhealthy plant indicates the presence of Ca. L. asiaticus infection. 2.The method of claim 1, wherein the miRNA399 sequence comprises TG X₁ X₂AAAGGAG X₃ X₄ TTGCC X₅ X₆ X₇ (SEQ ID NO:99), where X₁ is C or T, X₂ is Cor T, X₃ is A or C, X₄ is G, T, or A, X₅ is C or A, X₆ is T or G, and X₇is A or G.
 3. The method of claim 1, wherein the miRNA399 comprises asequence at least 90% identical to UGCCAAAGGAGAUUUGCCCGG (SEQ ID NO:9),UGCCAAAGGAGAGUUGCCCUA (SEQ ID NO:10), UGCCAAAGGAGAAUUGCCCUG (SEQ IDNO:11), and UGCCAAAGGAGAGUUGCCCUG (SEQ ID NO:12); the siRNA1005comprises a sequence at least 90% identical to ATAGATAATGGATCAACGGTTATA(SEQ ID NO:13); the siRNA1008 comprises a sequence at least 90%identical to TCGAACAAGGTAAGGATGTCA (SEQ ID NO:14) orCCTTGTTCGAACAAGGTAAGGATGTCATTCTTT (SEQ ID NO:100); and the siRNA1009comprises a sequence at least 90% identical to CTTCTAATAAACATGCATGAA(SEQ ID NO:15) or CGTCTTCTAATAAACATGCATGAACTTATT (SEQ ID NO:101).
 4. Themethod of claim 1, wherein the method further comprises detecting themRNA of a ubiquitin-conjugating enzyme E2 (UBC) gene.
 5. The method ofclaim 4, wherein the UBC mRNA comprises a sequence that is at least 90%identical to SEQ ID NO:84.
 6. The method of claim 1, wherein the methodfurther comprises measuring phosphate levels in the plant.
 7. The methodof claim 1, further comprising contacting the plant with phosphate or aphosphorus oxyanion if the plant is infected with Ca. L.
 8. A method oftreating HLB disease comprising contacting a plant infected with HLBdisease or having Ca. L. asiaticus-infection with phosphorus oxyanions,stimulating phosphate uptake in the plant by suppressing UBC mRNA or UBCpolypeptide and inducing phosphorus transporters in the plant, therebyameliorating symptoms of HLB disease in the plant.
 9. The method ofclaim 8, further comprising contacting the plant with phosphorusoxyanion solutions in sufficient amount to ameliorate symptoms of HLBdisease in the plant.
 10. A kit for detection of Huanglongbing (HLB)disease or Ca. L. asiaticus-infection, the kit comprising: one or moreagents that specifically detects miRNA399, siRNA1005, siRNA1008 orsiRNA1009.
 11. The kit of claim 10, further comprising oligonucleotideprimers comprising a sequence selected from the group consisting of SEQID NOs:85, 86, 87, 88, 89, 90, 91, 92, 93, 94, 95, and 96.